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I review the current status of the determination of helicity-dependent, or polarized, 
parton distribution functions from a comprehensive analysis of experimental data in 
perturbative quantum chromodynamics. I illustrate the latest achievements driven by 
new measurements in polarized proton-proton collisions at the Relativistic Heavy Ion 
Collider, namely the first evidence of a sizable polarized light sea quark asymmetry and 
of a positive polarized gluon distribution in the proton. I discuss which are the open 
issues in the determination of polarized distributions, and how these may be addressed 
in the future by ongoing, planned and proposed experimental programs. 
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Understanding how the nucleon spin is built up from the spin of quarks and gluons - 
and their orbital angular momentum - is one of the most challenging goals in hadron 
physicsP For instance, according to the Jaffe and Manohar helicity sum rul^ 

2 = -AE(/i 2 ) + AG(/r 2 ) + C q (n 2 ) + C g (v ?), (1) 

the one-half proton spin can be explicitly decomposed into contributions from quark 
and gluon spin, AE and AG, and from quark and gluon orbital angular momentum, 
C q and C g . Whether each of these terms allows for a unique field-theoretic definition 3 
in Quantum Chromodynamics (QCD) - and is possibly a gauge-invariant, physically 
meaningful, measurable quantity - has raised a major controversy in the last years, 
partly clarified only very recently.® 

Usually, Eq. 0 is probed by measuring spin asymmetries, i.e. differences of cross 
sections with opposite polarizations of initial-state particles, that arise in a large va¬ 
riety of high-energy, large-momentum-transfer processes.® Following factorization,® 
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these are described as a convolution between a short-distance part (that contains 
information on the hard interactions of partons in the form of process-dependent 
kernels) and a long-distance part (that contains information on the spin structure 
of the nucleon in the form of universal parton distributions). The former can be 
computed in perturbative QCD; the latter should be determined from experimental 
data; both depend on the factorization scheme and scale r. 

In this write-up, I review recent progress in the determination of helicity- 
dependent, or longitudinally polarized, Parton Distributions Functions (PDFs) 

A/(x,/x 2 ) = f{x,n 2 ) - f l (x,R 2 ), f = u,u,d,d,s,s,g, (2) 


defined as the momentum densities of partons with spin aligned along (^) or opposite 
(-1) the polarization direction of the parent nucleon. Polarized PDFs encode the spin 
structure of the nucleon, since these are related to the first two terms in Eq. Q: 


AS(^j 2 ) = [ dx [Aq(x, r 2 ) + Aq(x, r 2 )] A G(r 2 ) = f dxAg(x,R 2 ). 

q=u,d,s J ° J ° 


(3) 


The dependence of the PDFs on x, the momentum fraction of the proton carried 
by the parton, is genuinely nonperturbative, and, as such, must be inferred from 
data. These are usually supplemented with some theoretical constraints in order to 
achieve a meaningful PDF determination. First, PDFs must lead to positive cross 
sections: at leading order (LO), this implies that polarized PDFs are bounded by 
their unpolarized counterpart^] \Af(x,R 2 )\ < f(x,R 2 ). Second, polarized PDFs 
must be integrable: this corresponds to the assumption that the nucleon matrix 
element of the axial current for each flavor is finite. Third, it follows from SU(2) 
and SU(3) flavor symmetry that the first moments of the nonsinglet C-even PDF 
combinations, AT 3 = A u + — Ad + and ATg = Au + + A d + — 2As + (where A q + = 
Aq + Aq, q = u,d,s), are related to the baryon octet /3-decay constants, whose 
values are well measured P 

dxAT 3 = 1.2701 ± 0.0025 a s = [ dxAT a = 0.585 ± 0.025 . (4) 

Jo 

The dependence of the PDFs on r is perturbative. This obeys Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations^ and has been com¬ 
puted up to next-to-next-to-leading order (NNLO)j^jln principle, a determination 
of polarized PDFs at NNLO is then possible, even though based on inclusive DIS 
data alone: indeed, among the observables relevant for a determination of polarized 
PDFs, coefficient functions are known at NNLO only for the polarized DIS structure 
function rq jAl In practice, the impact of NNLO corrections is smaller than the cur¬ 
rent experimental uncertainties on data, hence moving from next-to-leading order 
(NLO) to NNLO determinations of polarized PDFs is not convenient so far. 



a Beyond LO, more complicate relations hold, but they have negligible impact on PDFsPl 
b In the polarized case, the computation has been completed only very recently ULJ 
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Several determinations of polarized PDFs of the proton (up to NLO and mostly 
in the MS factorization scheme) are presently available, aiming at unveiling how 
much large and uncertain are AS and AG, Eq. ([3]). Above all, they differ among 
each others for the procedure used to determine PDFs from data and for the included 
data sets (for details, see e.g. Chap. 3 in Ref. iTiJfl . Both are a source of uncertainty, 
but the former may be elusive: in this respect, the NNPDF collaboration has developed 
a methodology to reduce and keep it under control as much as possible. This is 
based on cutting-edge statistical tools, including Monte Carlo sampling for error 
propagation, neural networks for PDF parametrization, and closure tests for explicit 
characterization of procedural uncertainties.^ 

The bulk of the experimental information on polarized PDFs comes from neutral- 
current inclusive and semi-inclusive Deep-Inelastic Scattering (DIS and SIDIS) with 
charged lepton beams and nuclear targets. Because of the way the corresponding 
observables factorize, inclusive DIS data constrain the total quark combinations 
A< 7 + , while SIDIS data, with identified pions or kaons in the final state, constrain 
individual quark and antiquark flavors. In principle, both DIS and SIDIS data would 
constrain the gluon distribution A g via scaling violations, but in practice their effect 
is rather weak because of the small Q 2 range covered. 

Beside DIS and SIDIS fixed-target data, a significant amount of data from lon¬ 
gitudinally polarized proton-proton (pp) collisions at the Relativistic Heavy Ion 
Collider (RHIC) have become available recently^ though in a limited range of mo¬ 
mentum fractions, 0.05 < x < 0.4. On the one hand, longitudinal (parity-violating) 
single-spin and (parity-conserving) double-spin asymmetries for W ± boson produc¬ 
tion are sensitive to the flavor decomposition of polarized quark and antiquark 
distributions, because of the chiral nature of the weak interactions.^ On the other 
hand, double-spin asymmetries for jet and 7r° production are directly sensitive to 
the gluon polarization in the proton, because of the dominance of gluon-gluon and 
quark-gluon initiated subprocesses in the kinematic range accessed by RHICP^ 

Motivated by the interest in studying the effects of this piece of experimental in¬ 
formation, two new global analyses of polarized PDFs have been carried out in 2014, 
DSSV14 17 and NNPDFpoll. I. 18 These upgrade the corresponding previous analyses, 
DSSVOsfCS and NNPDFpoll . with data respectively on double-spin asymmetries 
for inclusive jet production 21 and 7r° production!^] and on double-spin asymmetries 
for high-pr inclusive jet productiorEUE^llMI and single-spin asymmetries for W± pro¬ 
duction. 25 Some new data by the COMPASS experiment have also been included 
in DSSV14 and NNPDFpoll. 1, respectively new DIS and SIDIS datdEH 27 and open- 
charrn leptoproduction dataP^ The new data have been included in NNPDFpoll. 1 
by means of Bayesian reweightingpl and in DSSV14 by means of a full refit. 

Overall, both DSSV14 and NNPDFpoll. 1 PDF determinations are state-of-the-art 
in the inclusion of the available experimental information. The data sets in the two 
analyses differ between each other only for fixed-target SIDIS and RHIC 7r° produc- 


Preliminary RHIC results included in Ref. 1191 have been replaced in Ref. 1171 with final results. 
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Figure 1. (Left) The polarized light sea quark asymmetry x(Au — Ad) from the NNPDFpoll.l 
and NNPDFpoll.1+ (supplemented with pseudodata for production at RHIC) PDF sets at 
Q 2 = 10 GeV 2 , compared to expectations from various models of nucleon structure!^ (Right) 
The polarized gluon xAg from the DSSV14 and NNPDFpoll. 1 PDF sets at Q 2 = 10 GeV 2 . 


tion measurements, included in DSSV14, but not in NNPDFpoll. 1. The information 
brought in by these data is complementary to that provided by RHIC W ± produc¬ 
tion and inclusive jet production data respectively, but this is less constraining.^ 
These data were not included in the NNPDFpoll. 1 analysis because fragmentation 
functions (FFs) enter the factorized expression of the corresponding observables: 
since FFs are nonperturbative objects on the same footing as PDFs, they are likely 
to introduce an additional source of bias in the PDF determination. The NNPDF 
methodology aims at reducing this bias as much as possible, hence the inclusion 
of these data would require the consistent determination of FFs within the NNPDF 
methodology, which is under consideration, though not yet available.^ 

The effect of RHIC data on the polarized PDFs of the proton is twofold. 

• The 2012 STAR data sets on W production,® 5 included in NNPDFpoll.l, pro¬ 
vide evidence of a positive A u distribution and a negative Ad distribution, with 
| Ad] > |Aii|. 18 The size of this flavor symmetry breaking for polarized sea quarks 
is quantified by the asymmetry Ad — Ad, which, in the NNPDFpoll.l analysis, 
turned out to be roughly as large as its unpolarized counterpart (in absolute 
value), though much more uncertain.^ Even within this uncertainty, polarized 
and unpolarized light sea quark asymmetries show opposite sings, with the polar¬ 
ized being definitely positive. This result starts to discriminate between different 
models of nucleon structure, see left panel of Fig. [T] specifically, some meson-cloud 
(MC) models are disfavored, while a more accurate experimental information is 
needed to establish whether chiral quark-soliton (CQS), Pauli-blocking (PB) or 
statistical (ST) models are preferred (these models are described in Ref. . 

• The 2009 STAR and PHENIX data sets on jet and 7r° production JHE2I included 
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in both DSSV14 and NNPDFpoll. 1, provide first evidence of a sizable, positive 
gluon polarization in the proton. A comparison of the gluon PDF in the two 
PDF sets is displayed in Fig. [I] (right panel). Comparable results, both central 
values and uncertainties, are found in the x region covered by RHIC data. The 
agreement between the two analyses is optimal in the range 0.08 < x < 0.2, where 
the dominant experimental information comes from jet data; a slightly smaller 
central value is found in the DSSV14 analysis, in comparison to the NNPDFpoll. 1, 
in the range 0.05 < x < 0.08, where the dominant experimental information 
comes from 7r° production data. Indeed, these are included in DSSV14 but are not 
in NNPDFpoll. 1. Nevertheless, best fits lie well within each other error bands, 
though NNPDF uncertainties tend to be larger than DSSV14 uncertainties outside 
the region covered by RHIC data. Very well compatible values of the integral of 
A 9, Eq. truncated over the interval 0.05 < x < 1, are found: at Q 2 = 10 

GeV 2 , this is 0.20±°;or for DSSV14JI3 and 0.23 ± 0.06 for NNPDFpoll. lP 

Despite the achievements described above, the lack of experimental data over a 
wide range of x and Q 2 values seriously limits the accuracy with which polarized 
PDFs can be determined. Several issues in our knowledge of the nucleon longitudinal 
spin structure are hence left completely open, as summarized below. 

• The size of the contribution of quarks, antiquarks and gluons to the nucleon spin, 
as quantified by their first moments, Eq. ([3]), is affected by large uncertainties. 
These come predominantly from the extrapolation into the small-a; region (x < 
10 -3 ), not covered by experimental data. In order to illustrate the situation, the 
running integrals of singlet and gluon distributions ( i.e . the quantities in Eq. ([3]) 
evaluated as a function of the lower limit of integration x m i n ) are displayed in 
Fig. [ 2 ] at Q 2 = 10 GeV 2 . It is apparent that, as x m i n decreases, the uncertainty 
on the integrals increases up to a size that prevents from any firm conclusion on 
their contribution to Eq. (JT]) . 

• A precision test of the Bjorken sum rule 3 ^ is presently not achievable within the 
accuracy of available data: indeed, a largely uncertain, and potentially substantial, 
contribution to it may arise in the small-a; region, x < 10 _3 P As a consequence, 
a determination of a s from the Bjorken sum rule is not competitive. 

• Fairly significant violations of SU(3) symmetry are advocated in the literature 
(see e.g. Ref. [34] for a review). In this case, an uncertainty on the octet axial 
charge, larger up to 30% than its experimental value, Eq. 0 , is found P It was 
shown 20 ] that a more conservative estimate of the uncertainty on as has a limited 
impact on the behavior of polarized PDFs in a global analysis. Nevertheless, if 
the octet sum rule in Eq. 0 were not imposed at all, data alone would not be 
sufficiently accurate either to constrain ATg or to test potential SU(3) violations. 

• Inclusive DIS data, together with nonsinglet axial couplings, Eq. Q, and kaon- 
tagged SIDIS data provide the sole constraint on the total strange distribution 
As + . No experimental information is available on individual As and As PDFs, 
and As = As is often assumed for phenomenological purposes. Mutually con- 
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Figure 2. The running integral f 1 dxAf(x,Q 2 ), f = E,</, from NNPDFpoll. at Q 2 = 10 

° ^min < < 

GeV 2 (upper panels), and their ratio to NNPDFpolEIC-E 3 ^based on EIC pseudodata (lower panels). 


sistent, sizable, negative values of the first moment of As + are found in both 
DSSlAtEHand NNPDfGHEH analyses, though they arise from rather different shapes 
of As+. This discrepancy may considerably depend on the kaon FF used to an¬ 
alyze SIDIS datapSl which are included in the DSSV determinations and are not 
in the NNPDF determinations. Further measurements of spin-dependent kaon pro¬ 
duction cross section in SIDIS and an improved determination of the kaon FF 
are needed to make any definitive conclusion. 

• Several models of nucleon structure have been developed, aiming at predicting 
the polarized PDF behavior at small and large x. In order to test whether data 
favor or unfavor them, computations of spin-dependent observables based on these 
models should be compared to the corresponding expectations based on a global 
analysis of polarized PDFs. Unfortunately, only a limited number of models are 
clearly incompatible with the predictions based on DSSV14 and NNPDFpoll. I. 3 ' 
More abundant and more accurate data at small and large x would be needed to 
discriminate among the other models. 

An intense experimental campaign is being devised to address these issues in the 
future. Some ongoing, planned and proposed prospects are discussed in the sequel, 
in order of increasing impact on the proton longitudinal spin structure. 

• The COMPASS experiment at CERN is about completing its experimental pro¬ 

gram dedicated to the longitudinal spin structure of the nucleon. It will provide 
additional high-precision results for the inclusive DIS structure function of the 
proton However, these are expected to be of limited impact in narrowing 

the small-a; extrapolation uncertainty on the full first moments, Eq. § , since 
they will extend only down to x ~ 0.004. 

• The CLAS and Hall-A experiments at JLAB have recently presented a large 
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amount of high-precision data respectively on the ratio of polarized to unpolarized 
proton and deuteron structure functions, j F\‘' d W^ and on the proton and 

neutron spin photoabsorption asymmetry, Al (’". 41 Since these data are taken in a 
kinematic region (large x , small Q 2 ) where dynamic higher-twist contributions to 
the Wilson expansion of g\ and resummation effects may be relevant J5HEU their 
inclusion in a global analysis of polarized PDFs is highly nontrivial, and not yet 
performed. The kinematic reach of the existing JLAB data is expected to extend 
to twice smaller x as well as to larger x values after the 12 GeV electron beam 
energy upgrade : 44 large luminosity and high resolution available will allow for a 
substantial reduction of PDF uncertainties in the medium-to-large x region. 

• The STAR and PHENIX experiments at RHIC are expected to provide high- 
impact results in the near term. l4 J As for W ± production, much smaller uncer¬ 
tainties on the single-spin asymmetry will be reached thanks to the substantially 
increased statistics from run-2013. Supplementing NNPDFpoll.l with W ± pro¬ 
duction pseudodata based on projected uncertainties will reduce the uncertainty 
on the polarized light sea quark asymmetry by a factor two, see Fig. [l] (the PDF 
determination including pseudodata is labeled NNPDFpoll. 1+). As for 7r° and jet 
production, future RHIC measurements up to run-2015 at center-of-mass energy 
yfs = 510 GeV will be sensitive to the polarized gluon PDF down to x ~ 10 -3 . 
Inclusion of n 0 and jet production pseudodata based on projected uncertainties 
in the DSSV14 analysis will reduce the uncertainty on A g by a factor twoP^l Com¬ 
plementary information on Ag will be provided by correlation measurements, e.g. 
in di-jet and di-hadron production processes. 

• A future high-energy, polarized Electron-Ion Collider (ElCjP^I will likely be the 
only facility to address all the above questions with the highest precision. The 
extension of the kinematic reach down to x ~ 10~ 4 and up to Q 2 = 10 4 GeV 2 
will allow for an accurate determination of A g via scaling violations in inclusive 
DIS, of Art and Ad via inclusive DIS at high Q 2 mediated by electroweak bosons, 
and of As via kaon-tagged SIDIS. Using simulated pseudodata at the eRHIC 
realization of an EIC,^ the impact of some of these measurements has been 
recently studiedES 4 ' It was found that the running integrals AE and AG will be 
determined with an accuracy of about respectively ±20% and ±10%, see Fig. [2] 
Would the EIC data confirm the DSSV14 or NNPDFpoll.l best fit behaviors of 
AE and AG, see Fig. [2j only a small fraction of the proton spin is expected to 
come from orbital angular momentum, C q ± C g in Eq. 0 ’ at Q 2 = 10 GeV 2 . 

A complete understanding of the origin of the proton spin is still lacking. Brand 
new facilities, such as the EIC, will be eventually required to elucidate how this 
is built up from the interplay between the intrinsic properties and interactions of 
quarks and gluons. Exploring the proton spin may lead to unexpected discoveries in 
the future, like it happened in the past for the evidence of nucleon substructure in 
the measurement of the proton anomalous magnetic moment; ultimately, this may 
shed light on confinement, an essential step towards our understanding of QCD. 
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